DNA unzipping and the unbinding of directed polymers in a random media. 
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We consider the unbinding of a directed polymer in a random media from a wall in d = 1 + 1 
dimensions and a simple one-dimensional model for DNA unzipping. Using the replica trick we show 
that the restricted partition functions of these problems are identical up to an overall normalization 
factor. Our finding gives an example of a generalization of the stochastic matrix form decomposition 
to disordered systems; a method which effectively allows to reduce dimensionality of the problem. 
The equivalence between the two problems, for example, allows us to derive the probability distri- 
bution for finding the directed polymer a distance z from the wall. We discuss implications of these 
results for the related Kardar-Parisi- Zhang equation and the asymmetric exclusion process. 



B 

i 

-a 
c 
o 

(N 
> 

in 
o 
\o 
o 

-I— > 
c3 



i 

C 

O 

o 



X 



The problem of a directed polymer in a random me- 
dia (DPRM) has received much attention for more than 
two decades 0. This stems from many reasons: It is a 
relevant model for a single vortex in a disordered type-II 
superconductor and it is one of the simplest examples 
of disordered systems for which, in some cases, exact re- 
sult can be obtained. Moreover, the model maps both 
to the Kardar-Parisi-Zhang (KPZ) equation which is 
perhaps the simplest non-linear stochastic growth equa- 
tion, and to the noisy Burger's equation jj]. The lat- 
ter, in one dimension, describes the long-time and large 
length-scale behavior of asymmetric exclusion processes 
(ASEP) 0, Q, which have been studied as prototypes 
of non-equilibrium systems 0, II] • The relations between 
the different models has been extremely fruitful: in some 
cases results which are hard to derive in one model can 
be easily obtained using another. 

Here we present a new, more subtle, relation between 
the DPRM problem and the DNA unzipping problem . 
In particular, we consider (a) the depinning of a DPRM 
from an attractive wall in d — 1 + 1 dimensions and (b) 
the force induced unzipping of a directed elastic line (in 
d > 1 + 1 dimensions) from a disordered columnar defect. 
We study a low temperature limit where the only excur- 
sion of the line from the defect occurs at the edge of the 
sample where the external force is acting (see Fig. 0. 
Both problems are disordered. In the DPRM problem 
the half-plane is taken to have uncorrelated point disor- 
der, while in the unzipping problem there is uncorrelated 
point disorder localized on the columnar defect. In the 
DPRM problem it is well known that as the strength 
of the disorder grows there is an unbinding transition 
whereby the polymer leaves the wall. Similarly, it is 
known that in the unzipping problem as the force act- 
ing on the line exceeds a critical value, the line leaves the 
potential. The unzipping model has been recently used 
in the context of single molecule experiments performed 
on DNA |9|, [lfj, [llj and also in relation to magnetic force 
microscopy experiments in type-II superconductors [T^ . 

Both problems have been treated in some detail before 
and superficially bear little resemblance. For example, 
the unbinding of the DPRM is a continuous phase tran- 
sition while the unzipping problem is a first-order phase 



transition. More interesting, the dimensionality of the 
problems is different. As we show the unzipping problem 
is effectively one-dimensional, while the DPRM problem 
is two dimensional. Nevertheless, we find that below the 
unbinding/unzipping transitions the replicated partition 
function of the DNA unzipping model is identical to the 
replicated generating function for the localization length 
of the polymer near the wall. Such an identity allows 
us to apply all known results from the unzipping prob- 
lem to the DPRM unbinding problem. For example, we 
can derive the distribution function of the DPRM local- 
ization length near the wall below the transition. This 
identity also allows one to perform large scale numerical 
simulations of the DPRM by effective reduction of the 
dimensionality of the problem. 
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FIG. 1: The two problems discussed in the text, (a) Unbind- 
ing of a line from an attractive wall in the presence of point 
disorder in d = 1 + 1. (b) Unzipping of a line with point 
disorder confined to an attractive line. 



We also discuss the implication of our result to both the 
KPZ equation and the ASEP. When the DPRM model 
is mapped to the KPZ equation, the wall acts as a free 
boundary condition for the growing interface and the at- 
tractive potential reduces the growth locally near the in- 
terface. The transition is approached by increasing the 
local growth rate at the interface. The relation to the 
unzipping problem allows us to derive the distribution 
function for the height profile near the transition. When 
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the DPRM is mapped to an ASEP the boundary corre- 
sponds to particles being injected at the end of a semi- 
infinite system with a rate which is slower than the par- 
ticle hopping rate in the bulk of the system. Here the 
transition is approached by increasing the injection rate 
of particles. Using our results we discuss the structure 
of the particle density near the transition. The relation 
between the models suggests a deeper relation which may 
not be transparent in our derivation. 

Equivalence between a two-dimensional (2d) unbinding 
problem and a one-dimensional (Id) unzipping problem 
is quite remarkable because lower dimensional problems 
are generally much easier to deal with. Without disorder 
our results can be understood in the context of stochastic 
matrix form (SMF) decomposition 0] as follows. First, 
we map the partition function of a 2d classical prob- 
lem (DPRM) into the partition function of a quantum 
Id problem, where the coordinate along the polymer is 
replaced by an imaginary time. Next we interpret the 
ground state wave- function of this quantum problem as a 
partition function of another classical Id problem (DNA 
unzipping). The latter step is called SMF decomposition. 
Here we present an example extending this method to a 
disordered problem which corresponds to a Hamiltonian 
which depends on (imaginary) time. We believe that our 
findings can be further generalized to a wider class of 
disordered problems. 

Unzipping of an elastic line from a columnar pin. Let 
us consider an external force f pulling the top end of 
a line from an attractive disordered columnar pin (see 
Fig. lb) . Neglecting excursions of the bound part of the 
line into the bulk, we can write the partition function as 
a sum over positions z where the line leaves the pin. The 
contribution of the unzipped part can be easily computed 
from the elastic energy [l2|: 



F {z)= -f -r(0) +[dz' 
Jo 



(1) 



Here r(z) denotes transverse coordinates of the line and z 
is the coordinate parallel to the defect. Integrating over 
all possible paths give s the corresponding free energy: 
Fq(z) = — / 2 z/27 19]. The free energy of the bound 
part, up to aconstant, is F\{z) — Vo(L — z)+ ff dz'U(z'), 
where Vb is the value of the attractive potential, L is the 
length of the columnar defect which is assumed to be 
very large, and U (z) is a random uncorrelated potential 
with zero mean satisfying U(zi)U(z2) — AS(zi — z-i). 

The partition function Z is a sum of the corresponding 
Boltzmann weights over all possible values of z 
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Z(z) 
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For simplicity we work in the units, where ksT = 1. In 
the equation above e = (/| — f 2 ) /2j, where / is the force 
applied to the end of the flux line and f c = Vb| 
is the critical force. Note that the partition function 
@ describes effectively a one-dimensional problem. The 
statistical properties of the unzipping transition can be 
obtained by replicating Eq. @ [lj] : 



Z n = Z^= dz x ... dz n c-^=i F ^\ 
Jo Jo 
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where the overline denotes averaging over point disorder. 
The averaging procedure can be easily done for a positive 
integer n. First we order the coordinates Zj, where the j th 
replica unbinds from the pin according to: < z\ < z 2 < 
• • • < z n . Then for z G [0, z\) there are no replicas bound 
to the columnar pin, for z G [z\, z-j) there is one replica 
on the pin until finally for L > z > z n all n replicas are 
bound to the pin. Using this observation and explicitly 
averaging over the point disorder in Eq. Q we arrive at: 



Z„ = n! 



- £ tz j -A/2j 2 (z j + 1 - Z] ) 

dz n . .. dz n e J - 1 (5) 



0<zi~-<z n <L 



where we use the convention z n+ \ = L. The integral 
above is straightforward to evaluate in the L — > oo limit: 
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_ /A/2L ( _2_V' r(l + l/n n - n) = etl 2 A/2i 



(6) 



where e„ = e + An and K n = A/2e„. The exponential 
prefactor is the contribution of the whole pin while the 
rest of the expression is the (L independent) contribution 
from the unzipped region. The disorder averaged free 
energy is given by the limit F — — lim n _ > o('^n — l)/ n - 
Using Eq. © one obtains 



(7) 



where ^(x) is the digamma function and k — A/2e. The 
unzipping transition occurs at e = or equivalently at 
k — > oo . The expression (JJJ is identical to the one found 
in Ref. |15| using the Fokker-Planck equation approach, 
supporting the validity of the replica calculation. 

Unbinding from a wall. Next we consider the unbind- 
ing of a directed polymer (or interface) from an attractive 
wall, in the presence of point disorder in the bulk of the 
system, in d = 1 + 1 dimensions (see Fig. [TJi) . The elastic 
free-energy of the problem is given by 



where up to an unimportant constant additive term VqL: 



F m {z) = F (z) + F x {z) = tz+ f dz'U(z'). 

J 2 



(3) 



F ub = / dr 



l(d T z(ryf + V(z(r))+^z(r),r) . ,M 



Here z(t) denotes the distance of the polymer from the 
wall at position t, 7 is the line tension, V(z) is a short 
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range attractive potential near the wall placed at z — 
and /z(z, r) is the contribution from the point disorder. 
The free energy of this problem was obtained first, using 
a replica calculation, by Kardar ^(|. For completeness, 
here we outline the main points of the derivation. The 
overall weight of paths connecting points (0, 0) and (z, r), 
W(z,t), can be calculated from 



d T W(z,T) 



[^(z,T) + V(z)- 1 d 2 z ] W(z,t) 
H(z,t)W(z,t). (9) 



After replicating the Hamiltonian n times one obtains 



H r , 



■E 

n = l 



V(z a )]-a S(zo 

a<j3 
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where we have assumed a Gaussian distribution of fi with 
zero mean and variance a. Averaging over disorder elim- 
inated the z dependence of /i. The ground-state wave 
function ip (as well as the energy) of the problem were 
obtained by Kardar using the Bethe ansatz. We skip the 
details as they can be found in Ref. . For the permuta- 
tion P of particles such that < zp\ < zp2 < . . . < zp n 
and below the transition ip ~ exp (— ^™ =1 XaZp a )- Here 
Xa = A + 2(q — l)x, X — cr/47 and A depends on the 
strength of the attractive potential at the wall. We note 
that the wave function gives the probability distribution 
of the directed polymer at the top of the sample (r = L), 
and thus it is proportional to the restricted partition 
function W(L, z). Integrating ip over all coordinates gives 
the generating function of moments of z |Tfjj : 



N n = / dz 1 



dz n ip = 



1 r(l + l/« n -n) 

x n r(i + i/ K „) : 



(ii) 



where now K n = x/ e n an d £n = A — x + ^X n - Note that 
with a proper identification of e n and n n 



= N„ 



(12) 



This is the main result of our paper. In fact, one can 
check that not only the generating functions Q n and N n 
coincide, but also the wave function of the DPRM prob- 
lem is identical (up to an unimportant constant) to the 
restricted partition function of the DNA unzipping prob- 
lem. This equivalence implies that below the transition all 
the moments and cumulants of the unzipped length and 
of the distance of the polymer from the wall are equiv- 
alent. In particular, the weight of finding the DPRM a 
distance z from the wall for a particular realization of dis- 
order is given by Z(z), defined in Eq. @. Note that the 
free-energies of the two models are distinct. Indeed Q n 
contains information on the free-energy of the unzipping 
problem. However, for the DPRM N n only describes the 
spatial distribution of z at the top of the sample. In par- 
ticular, the sum over weights of a DPRM, which ends at 
the point z at the top of the sample, W(z, L), is given by 



W(z,L) = e 



V>(*) 



Jdz'ip{z'] 



= e~ E ° L 



Z(z) 



(13) 



where Z(z) is defined in Eq. JSJ, Eq is a sample de- 
pendent free-energy and Z is the normalization factor. 
We note that the equivalence between ip{z) an( i Z( z ) 
gives an example of the SMF decomposition 0] for 
a disordered DPRM problem corresponding to a time- 
dependent Hamiltonian. 

It is easy to see that for both models 



(z) =dF/de={Ke)- 1 ¥ 1 \l/n), 



(14) 



where fW(i) stands for the n-th derivative of the 
digamma function. The expression above predicts a 
crossover from (z) w 1/e for k *C 1 (far from the transi- 
tion) to (z) w k/c = A/e 2 for k ^S> 1 (close to the transi- 
tion) as was noted previously for both unzipping .9] and 
unbinding 16] problems separately. Similarly 



w = (z 2 ) - (z) 2 = d 2 F/de 2 = -{k€)- 2 ¥ 2 \\/k). (15) 

Here there is a crossover from w ~ 1/e 2 for k <C 1 to 
w ~ 2k/ e 2 = A/e 3 for k» 1. As has been noted in the 
context of DNA unzipping y / w/(z) changes from being 
unity for k -C 1 to ~ e 1 / 2 for k ^> 1. Thus close to the 
transition, thermal fluctuations become negligible. 

Calculation of the second moment. With a little more 
work we can evaluate higher moments of the distribu- 
tion. In particular, the second moment, which gives the 
variance. To do this we consider the generating function: 



W n =n\ dzi 



dz r , 



- E e j z j +A/2j 2 (z j - 1 -z j ) 



(16) 



0<2 Tl ■ ■ ■< z\ <oo 



The second (and similarly the higher) moments can be 
found by differentiating W„ with respect to ef. 



1 1 



n 0,1 



(- 2 > = ]i-— :E 



d 2 w n 



oWnnfj de 2 

3 = 1 ■? 



(17) 



Evaluating the simple integral in Eq. (|16fl we find that 



w„ = i/n;u ELi^- aj 2 /2 



and correspondingly 



(z 2 ) 



1 1 n 2 ™ 1 
= — lim — y y — . 

e 2 n^o n 1 — ki ^— ' Ml — kk) 



(18) 



This sum can be evaluated using a trick similar to the 
one suggested by Kardar fl6l |: 



2k 2 



dxdy 



1 



ye" 



K v 4. e 2 y e Kx ~ x 



- [e«v + 



x>y>0 



(C + *(!/«)) ; 



(19) 



where C « 0.577 is the Euler's constant. In the limit of 
weak disorde r or high temperature n <C 1, not surpris- 
ingly, we get (z 2 ) pa 2/e 2 , which agrees with the Poisso- 
nian statistics of z with a average given by (z) — 1/e. In 
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the opposite limit k 3> 1 one finds (z 2 ) = An 2 /e 2 — 4(z) . 
Thus in this limit the relative width of the distribution 
5z/(z), where Sz 2 — (z 2 ) — (z) 2 , is larger by a fac- 
tor of -\/3 than that in the high temperature regime. 
The distribution becomes superpoissonian at large k. 
In fact at k — ► oo one can derive the full distribution 
function P K _ >00 (z) using extreme value statistics pi IT^: 
V^ooiz) w e/nG(ze/n), with 

G(x) = -L e"*/ 4 - ierfc(Vi/2), (20) 

where erfc(x) is the complimentary error function. It is 
easy to check that this distribution indeed reproduces 
correct expressions for the mean and the variance. We 
additionally performed direct numerical simulations of 
the partition function and got excellent agreement 
with predictions of Eqs. JTSJl and 0- 

The quantity (z 2 ) is not always experimentally acces- 
sible. For example, in the unzipping experiments it is 
easier to measure thermal average, (z), in each experi- 
mental run. Then the variance of the distribution will 
be characterized by (z) 2 . The difference between the two 
expectation values is given by w found in Eq. I|15|) . Defin- 
ing b~z\ = (z) 2 — (z) and using Eqs. (|19(l and i|15fl we 
find that 5zt/(z) ~ \/k/2 in the weak disorder limit 
(k < 1) and Sz T /(z) sa \/3 - 1/(%/3k) for k > 1. We 
plot both Szt and Sz versus the disorder parameter n 
in Fig. [3 Note that importance of the order of thermal 




FIG. 2: Dependence of the relative width of the distribu- 
tion on the disorder parameter k. The two curves correspond 
to different averaging over temperature and disorder (see the 
text for details). The horizontal line at v3 denotes the asymp- 
totic value of both 5z/{z) and Szt/{z) at k — > oo 

and disorder averaging also appears in the calculation 
of higher moments of z becoming irrelevant only in the 
zero-temperature (strong disorder) limit (k — > oo). 

Implications. We now turn to discuss the implications 
of our results for the KPZ equation and the ASEP. The 
height variable in the KPZ equation is well known to 



be related to the DPRM through the Cole-Hopf trans- 
formation h(z,r) = — ln(W(z, r)). The latter together 
with Eq. @ yields 

d T h(z, t) = -fd 2 z h(z : t) - l{d z h(z, t)) 2 + V(z) + n(z, r) . 

Now r represents a time coordinate for the growing in- 
terface. The pinning potential at the origin leads to a 
reduced growth rate at a free boundary of the interface 
(for more details on the correspondence see, for exam- 
ple, As stated previously, the relation between the 
unzipping problem and the DPRM does not give the full 
information about W(z,t). It only contains its spatial 
behavior but misses the prefactor (see Eq. H13[) h which 
corresponds in the KPZ picture to the overall height of 
the interface. However, the mapping implies that up to 
this overall height the interface is described by Eq. 10, 
namely h(z,r) — h a (r)+F uz (z), where h a {r) is the mean 
height and F uz (z) is a tilted random walk. The proba- 
bility for an interface profile F uz (z) is equal to the prob- 
ability of drawing the random contributions U(z) (see 
Eq. ©). The correspondence implies, for example, that 
h{z) — h(0) has a Gaussian distribution with a variance 
which increases linearly with z. The localization of the 
polymer near the unbinding transition corresponds to a 
single dominant minima in the free energy (equivalent 
to a minimum in the interface), the average distance of 
these minima from the interface behaves as 1/e 2 (see the 
discussion after Eq. ((Hfl h 

Next, we turn to implications of our results for the 
ASEP. We consider a semi-infinite system where parti- 
cles hop into a one dimensional lattice at the left end 
with a rate a and hop in the bulk of the system to the 
nearest neighbor on the right with a rate 1. In the con- 
tinuum limit the behavior of the system is captured by 
the equation 

d rP {z) = -^d 2 z p(z) + 2p(z)d z p(z)-d z p(z) + U(x)+T 1 (z,T) . 

(22) 

where U (x) represents injection of particles into the sys- 
tem and rj{z, r) = d z p(z, r) (with p(z, r) as before drawn 
from a Gaussian distribution) is a conserving noise. The 
equation can be obtained from Eq. I|21|l by setting 7 = 
1/2 and defining a variable m(z, t) = (z—h(z 7 T))/2 (note 
that d z h(z,r) satisfies the noisy Burgers equation). It 
is straightforward to verify that p(z,r) — d z m(z,T) in- 
deed satisfies Eq. (|2"2"|) . Using the same reasoning as in 
the KPZ equation we find that in the steady state the 
density profile is described by p(z) = (1 — e — U{z))/2. 
This result suggests that in a discrete realization of the 
model a particle appears at a given site with a constant 
probability and there are no correlations between differ- 
ent sites. Our results hold, in ceratin limits, even when 
another boundary is included in the system. For exam- 
ple when at the other end the particles are ejected at a 
high rate. Indeed it is known that for the ASEP with 
open boundaries the density is flat near the end where 
particles are injected into the system 0,0. However, we 
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emphasize the lack of correlations near the left boundary 
arbitrarily close to a continuous transition. Similarly, for 
the interface growth with two boundaries our results ap- 
ply if the growth rate is sufficiently slow only near one 
of the walls. We have performed numerical simulation 
which verified these statements. These results will be 
presented elsewhere [l8|. 

In conclusion, we demonstrated that there is exact 
mapping between the partition function of the DNA un- 
zipping transition and the spatial distribution of a DPRM 
unbinding from a wall. This mapping allowed us to apply 
some known and newly derived results of the simpler un- 



zipping problem to the DPRM problem. We also showed 
how this mapping can be used to derive results for the 
KPZ equation near a boundary and about asymmetric 
exclusion process with open boundaries. 
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